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Abstract: We give a probabilistic introduction to determinantal and per- 
manental point processes. Determinantal processes arise in physics (fermions, 
eigenvalues of random matrices) and in combinatorics (nonintersecting paths, 
random spanning trees). They have the striking property that the number 
of points in a region D is a sum of independent Bernoulli random variables, 
with parameters which are eigenvalues of the relevant operator on L'^{D). 
Moreover, any determinantal process can be represented as a mixture of 
determinantal projection processes. We give a simple explanation for these 
known facts, and establish analogous representations for permanental pro¬ 
cesses, with geometric variables replacing the Bernoulli variables. These 
representations lead to simple proofs of existence criteria and central limit 
theorems, and unify known results on the distribution of absolute values in 
certain processes with radially symmetric distributions. 
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Fig 1. Samples of translation invariant point processes in the plane: Poisson (left), determi- 
nantal (center) and permanental for K{z,w) = 2 +1^1 ), Determinantal processes 

exhibit repulsion, while permanental processes exhibit clumping. 


1. Introduction 

Determinantal point processes were first studied in 1975 by Macchi who was 
motivated by fermions in quantum mechanics. Determinantal processes arise 
naturally in several other settings, including eigenvalues of random matrices, 
random spanning trees and nonintersecting paths; see, e.g., Burton and Peman- 
tle 1^, Soshnikov EH. Lyons EHI, Lyons and Steif ED. Shirai and Takahashi Ea, 
Johansson El. Borodin, Okounkov and Olshanski |2|, and Diaconis [S]. A de¬ 
terminantal point process, on a Polish space A (assumed locally compact) with 
a reference measure (assumed Radon) n, is determined by a kernel K(x,y): the 
joint intensities of the process can be written as det{K(xi,Xj)). The kernel de¬ 
fines an integral operator 1C acting on L^(A) that is assumed to be self-adjoint, 
non-negative and locally trace class, i.e. for every compact D the eigenvalues, 
{Af}, of the operator 1C restricted to D satisfy < oo. Determinantal 

point processes have a special property (Shirai and Takahashi m Proposition 
2.8) that seems to have only been used in special cases (P, [^1: 

In a determinantal process, the number of points that fall in a compact set 

D C A, has the same distribution as a sum of independent Bernoulli(\f') 

random variables. 
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The proof is immediate from well known formulas for the generating function 
of particle counts. However, we give a proof starting from hrst principles and 
avoid the use of Fredholm determinants. 

Permanental processes are defined analogously and are the counterpart of 
determinantal processes for modeling bosons (see Macchi | 22 )- They fall into 
the more general class known as Cox processes |n]- In this case we have: 

In a permanental process, the number of points that fall in a compact set 
D C A, has the same distribution as a sum of independent geometric() 

random variables. 


For the examples of interest to us, operator theoretic nuances are not essen¬ 
tial. The reader will not miss anything significant by keeping in mind just the 
following cases. 

1. A is a finite set and AT is a Hermitian non-negative definite |A| x |A| matrix, 
and fi is the counting measure on A. 

2. A is an open set in R‘^, pL is Lebesgue measure and K{x,y) is continuous 
function defining a self-adjoint non-negative integral operator K. on if {A). 

A point process in a locally compact Polish space A is a random integer¬ 
valued positive Radon measure X on A. (Recall that a Radon measure is a 
Borel measure which is finite on compact sets.) If X almost surely assigns at 
most measure I to singletons, it is a simple point process; in this case X can be 
identified with a random discrete subset of A, and X{D) represents the number 
of points of this set that fall in D. 

The distribution of a point process can, in most cases, be described by its 
joint intensities (also known as correlation functions). 

Definition 1. The joint intensities of a point process X w.r.t. p are functions 
(if any exist) pk '■ A^ ^ [0, oo) for A: > 1, such that for any family of mutually 
disjoint subsets Di,... ,Dk of A, 


E 








Pk{xi,Xk)dp{xi) ... dp{xk), 


( 1 ) 


Remark 2. For overlapping sets, the situation is more complicated. Restricting 
attention to simple point processes, pk is not the intensity measure of X^, but 
that of X'^^, the set of ordered A:-tuples of distinct points of X. Indeed, CJ 
implies (see HHUHEg) that for any Borel set R C A^ we have 

E#(R n A^'") = / pkixi,...,Xk)dp{xi)...dp{xk). (2) 

J B 


When R = n for a mutually disjoint family of subsets Ri,..., of A, 
and k = X]i=i I^e left hand side becomes 
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For a general point process X, observe that it can be identified with a simple 
point process X* on A x {1, 2, 3,...} such that X*{D x {1, 2, 3,...}) = X{D) for 
Borel D C A. Thus, if X{D) has exponential tails for all compact D C A, then 
the joint intensities determine the law of see mmsi. In particular, Theorems 
13 and EH below imply that this is the case for determinantal and permanental 
processes governed by a trace class operator, since a convergent sum of Bernoulli 
(or geometric) variables always has exponential tails. 

Definition 3. A point process A on A is said to be a determinantal process 
with kernel K if it is simple and its joint intensities satisfy: 


Pk{xi,...,Xk) = det{K{xi,Xj))^^^ -^i^, (4) 

for every k > 1 and xi,... ,Xk € A. 

Remark 4. We postulate determinantal processes to be simple because we have 
adopted equation o as the definition of joint intensities. If instead, we start with 
the slightly more restrictive definition of joint intensities as explained above in 
RemarkEl it follows that the process must be simple. Nevertheless, postulating 
simplicity is more in tune with the conventions of Physicists who often consider 
a determinantal process with k points, not as a random counting measure but 
as a random point in A^/Diagonals, where “Diagonals” denotes the subset of 
points of A^ with at least two co-ordinates equal. This viewpoint together with 
a postulate on the behaviour of quantum amplitudes under exchange of particles 
leads naturally to Determinantal and Permanental processes (and additionally 
to “fractional statistics”, when A is two dimensional, see for instance M)- 

Consider a kernel K that defines a self-adjoint trace-class operator JC. Mac- 
chi m and Soshnikov m showed that there exists a determinantal point pro¬ 
cess with kernel K if and only if all eigenvalues of /C are in the interval [0,1]. In 
Section El we shall give a probabilistic proof of this fact. 

Remark 5. When we speak of a kernel AT on A^, a priori it is only defined only 
almost everywhere w.r.t. p, x p, and thus quantities like f K{x,x)dp{x) that 
appear in the definition of joint intensities are not defined. This can be made 
sense of as follows. 

The kernel K defines a self-adjoint integral operator 1C that has eigenfunctions 
(j)k and eigenvalues Afc, where (pk are orthogonal in A^(/r). In particular, there is a 
set Ai C A such that pk are all defined point-wise on Ai and /r(AJ) = 0. At least 
in the case when JC has finite rank, this shows that K(x,y) = J2k^k'Pkix)4)j.{y) 
is well defined on Ai x Ai, and that is sufficient to define / K{x, x)dp{x) etc. For 
more details on this point, see Lemmas 1,2 of the survey paper of Soshnikov m 

Recall that the permanent of an n x n matrix M is defined as 

n 

per(M) = ^ 

-K^Sn i—1 


( 5 ) 
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Definition 6. A point process A on A is said to be a permanental process 
with kernel K if its joint intensities satisfy: 

Pk{xi,...,Xk) = per {K{xi,Xj))j^^^.^i^, (6) 

for every k > 1 and xi,... ,Xk G A. 

For any kernel K that defines a self-adjoint non-negative definite operator 
K, (i.e., the eigenvalues of K, are nonnegative), there exists a permanental point 
process with kernel K. We shall give a proof of this known fact in Section 0] 
Now we state the main theorems. The most common application of the fol¬ 
lowing theorem is to describe the behavior of a determinantal process already 
restricted to a subset. 

Theorem 7. Suppose X is a determinantal process with a trace-class kernel K. 
Write 

n 

K (x, p) = J2 ’ 

k=l 

where 4>k o-xe normalized eigenfunctions of K with eigenvalues Xk G [0,1]. (Here 
n = oo is allowed). Let Ik, 1 < k < n be independent random variables with 
Ik ~ Bernoulli{Xk)■ Set 


Ki{x, y)=Y^ Ik(t)k{x)(j)k{y)- 

k=l 

Kj is a random analogue of the kernel K. Let Xj be the determinantal pro¬ 
cess with kernel Kj (i.e., first choose the Ip’s and then independently sample a 
discrete set that is determinantal with kernel Kj). Then 

X = Xi. (7) 

In particular, the total number of points in the process X has the distribution of 
a sum of independent Bernoulli(Xk) random variables. 

In the special case of random spanning trees of a finite graph, Bapat ^ 
was the first to observe the last fact stated above. Namely, he proved that the 
number of edges of the spanning tree falling in a subset of edges of the given 
graph has the distribution of a sum of independent Bernoullis. 

When the Dfs are related in a special way, there exists a simple probabilistic 
description of the joint distributions of the counts X{Di). 

Definition 8. Given an integral kernel K acting on lf{K), the subsets Di,..., Dr 
of A with D = UiDi are simultaneously observable if the eigenfunctions of 
the restricted kernel K\dxd are also eigenfunctions of Kl^ixOi acting on L‘^{Di) 
for every 1 < i < r. 

The motivation for this terminology comes from quantum mechanics, where 
two physical quantities can be simultaneously measured if the corresponding 
operators commute. 
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Proposition 9. Under the assumptions of Theorem^ let Di C A, 1 < i < r 
be mutually disjoint and simultaneously observable. Let be the standard basis 
vectors in K.’’. Denote by 4>k, the common eigenfunctions of K on the Di’s and 
by Xk^i the corresponding eigenvalues. Write Xk = ^k,i and note that Afc < 1. 
Then 

(A’pi), . . . , X{Dr)) = . . , 6,r) , (8) 


where = (Cfe.ii ■ ■ • Ak,r) are independent for different values of k, with P(^fe = 
Gi) = Xk,i for 1 < i < r and P(|fe = 0) = 1 — Xk- In words, {X{Di),... ,X{Dj.)) 
has the same distribution as the vector of counts in r cells, if we pick n balls and 
assign the k^^ ball to the i^^ cell with probability Xk,i (there may be a positive 
probability of not assigning it to any of the cells). 

Theorem 10. Suppose X is a permanental process in A with a trace-class kernel 
K. Write K(x,y) = YIk=i^k4>k{x)4>i,(]j), where (fk are normalized eigenfunc¬ 
tions of K with eigenvalues Xk (n = oo is allowed). Let a = (oi,..., a„), where 
ai are non-negative integers such that i = £{a) = oi + • • • + < oo and let Z°‘ 

be the random vector in with density: 


p{zi,...,ze) = 


1 


£!ai! • • • Ora! 


per 


{<>^’1(2:1) 


</)i(z£)}ai 


( 9 ) 


where the notation {fi^zf) ■ ■ ■ <fii{zi)'\ ai indicates that the row 4>i{zi) ■ ■ ■ 4’i{zt) is 
repeated at times. Let X^ be the point process determined by Z", i.e., Xs{D) is 
the number of j < £ such that € D. Let 71 ,... , 7 ^ be independent geometric 

random variables with V(yk = s) = ^ > for s = 0,1, 2,.... Then 

Zv. — ^ 

where 7 = ( 71 ,..., 7 „). In particular, X{A) has the distribution of a sum of 
independent geometric{j^^) random variables. 

Remark 11. The density given in m has physical significance. Interpreting 
the functions (fk as eigenstates of a one-particle Hamiltonian, 0 gives the dis¬ 
tribution for £ non-interacting bosons in a common potential given that ai of 
them lie in the eigenstate fi. This density is the exact analogue of the density 


p{zi,...,zi) = 


£>. 


det 






( 10 ) 


which appears in Theorem [3 and gives the distribution for £ non-interacting 
fermions in a common potential given that one fermion lies in each of the eigen¬ 
states ,..., ((> 7 . The fact that 03 vanishes if a row is repeated illustrates 
Pauli’s exclusion principle, which states that multiple fermions cannot occupy 
the same eigenstate. See |H] for more details. 
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Theorem 12, Under the assumptions of Theorem [23 suppose Di,..., Dr are 
simultaneously observable as in definition\^ Denote by (j)i~ the common eigen¬ 
functions of K on the Di’s and by Xk^i the corresponding eigenvalues. Then 

{XiDi), X{Dr)) = ^ (77fc,l, . . . , Pk,r) , (11) 

k 

where {rjkp, ■ ■ ■ ,r]k,r) are independent for different values of k, for each k, the 
sum rjk = Vk,i has a geometric distribution with mean Xk '■= and 

given Vk,i = N, 

{rjk 1 , - ■ ■ ,Vkr) = Multinomial f N; . ■ 

V 



2. Determinantal processes 


We begin with a few important examples of determinantal processes. 

Example 13 (Non-intersecting random walks). Consider n independent 
simple symmetric random walks on Z started from ii < i 2 <■■■< in where 
all the ij’s are even. Let Pi jit) be the t-step transition probabilities. Karlin 
and McGregor lEl show that the probability that the random walks are at 
ji < 32 < ■ ■ ■ < jn at time t and have mutually disjoint paths is 


det 


r^jiit) ... Piljni^) 
r„jiit) ... Pinjnit) 


It follows easily that if t is even and we also condition the walks to return to 
ii,..., in at time t, then the positions of the walks at time t/2 are determinantal 
with a Hermitian kernel. See Johansson M for this and more general results. 

Example 14 (Uniform spanning trees). Let G be a finite undirected graph 
and let E be the set of oriented edges (each undirected edge appears in E with 
both orientations). Let T be uniformly chosen from the set of spanning trees 
of G. For each directed edge e = vw, let X® := — '^wv denote the unit flow 

along e. Define 

H = {/ : E ^ R : f{vw) = —f{wv) Vu, w} 

★ = spanj^X'"'" : where u is a vertex.} 

W 

n 

<0 = span{^^ X®’ : ei,..., e„ is an oriented cycle}. 

i=l 


It is easy to see that H = ★©<(>. Now, define P := the orthogonal 

projection onto ★. Then, the set of edges in T forms a determinantal process 
with kernel K{e,f) := with respect to counting measure on the set of 
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edges. This was proved by Burton and Pemantle who represented K{e,f) 
as the current flowing through / when a unit of current is sent from the tail to 
the head of e. The Hilbert space formulation above is from BLPS [2]. 

Example 15 (Ginibre ensemble). Let Q be an n x n matrix with i.i.d. 
standard complex normal entries. Ginibre m proved that the eigenvalues of Q 
form a determinantal process in C with the kernel 


Kn{z,w) 


ie-i(hl'+l“l') 

TT 


n—1 


E 


(zW)'^ 

k\ 


with respect to Lebesgue measure. As n ^ oo, we get a determinantal process 
in the plane with the kernel 


K{z,w) = —e 

TT 


1 -i(P|2+i«,i= 


oo / _s u 

,^{zwY 


fe =0 


k\ 




Example 16 (Zero set of a Gaussian analytic function). The power series 
{i{z) where a„ are i.i.d. standard complex normals, defines a 

random analytic function in the unit disk (almost surely). Peres and Virag m 
show that the zero set of fi is a determinantal process in the disk with the 
Bergman kernel 


K{z, w) 


1 

7r(l — zw)^ 


1 

-Y^{k + l){zw)^ 

TT ^' 

k=0 


with respect to Lebesgue measure in the unit disk. 


Determinantal projection processes: motivation and construction 


The most general determinantal processes are mixtures of determinantal projec¬ 
tion processes, i.e. processes whose kernel Kh defines a projection operator JCh 
to a subspace H C L^(A) or, equivalently, KH(x,y) = X]fc=i 4>k(,x)4>ki.y) where 
{(j)k} is any orthonormal basis for H. 


Lemma 17. Suppose X is a determinantal projection process on A, with kernel 
K{x,y) = 'Y^^=i4>k{x)4)^.{y) where {(jk} is an orthonormal set in LS{K). Then 
the number of points in X is equal to n, almost surely. 


Proof. The conditions imply that the matrix {K{xi, has rank at most 

n for any fc > 1. From Q, we see that E = 0 for /c > n. This shows 
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that < n almost surely. However, 


E[A’(A)] 



A 


A 


pi{x) dfi{x) 


K{x, x) dfi{x) 



n. 


Therefore A (A) = n, almost surely. 


□ 


Despite the fact that determinantal processes arise naturally and many im¬ 
portant statistics can be computed, the standard Definition|21is lacking in direct 
probabilistic intuition. Below we present an algorithm that is somewhat more 
natural from a probabilist’s point of view, and can also be used for modeling 
determinantal processes. 

In the discrete case, the projection operator JCh can be applied to the delta 
function at a point, and we have K.h5x{‘) = K{-,x). In the general case we take 
this as the definition of IChSx- Let || • |j denote the norm of L^(/r).The intensity 
measure of the process is given by 


Ph{x) = pi{x)dp.{x) = \\JCH5x\\^dpL{x). 


( 12 ) 


When p is supported on countably many points, we have ||/Ci/i5a;|| = dist((5a;, H^) 
(where T denotes orthocomplement), giving a natural interpretation of the in¬ 
tensity Pi in general. 

Note that ph[^) = so pn/ dim{H) is a probability measure on A. 

We construct the determinantal process as follows. Start with n = and 

Hn = H. 

Algorithm 18. 

• If n = 0, stop. 

• Pick a random point Xn from the probability measure pH„/n. 

• Let Hn-i C Hn be the orthocomplement of the function JCn^dx in In 
the discrete case, iL„_i = {/ S Hn : f{Xn) = 0}. Note that dim(iL„_i) = 
n — 1 a.s. 

• Decrease n by 1 and iterate. 

Proposition 19. The points {Xi ,..., A„) constructed by A Inorithm X 1 fA are dis¬ 
tributed as a uniform random ordering of the points in a determinantal process 
X with kernel K. 

Proof. Let = TChSx, ■ Projecting to Hj is equivalent to first projecting to H 
and then to Hj, and it is easy to check that lCH,dxj = K.Hj'fj- Thus, by ifT^ . the 
density of the random vector {Xi, ..., A„) constructed by the algorithm equals 



J 
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Note that Hj = H Ci ..., and therefore V = 0^=1 is 

exactly the repeated “base times height” formula for the volume of the par¬ 
allelepiped determined by the vectors 'ipi,... ,'ipn in the finite-dimensional vec¬ 
tor space H C L^(A). It is well-known that equals the determinant of the 
Gram matrix whose i,j entry is given by the scalar product of V'iiV'j) is 
/ipiijjjdn = K{xi,Xj). We get 


p{xi,...,Xn) = ■^det{K{x^,Xj)), 

so the random variables Xi,..., are exchangeable. Viewed as a point process, 
the n-point joint intensity of {Xj}^^^ is nlp{xi,..., Xn), which agrees with that 
of the determinantal process X. The claim now follows by Lemma El D 

Example 20 (Uniform spanning trees). We continue the discussion of Ex¬ 
ample El LU Gn+i be an undirected graph on n -|- 1 vertices. For every edge e, 
the effective resistance when a unit of current is sent along e is R{e) = (I®, /®). 
To use our algorithm to choose a uniform spanning tree, proceed as follows: 

• If n = 0, stop. 

• Take Xn to be a random edge, chosen so that P(V„ = e^) = ^R{ei). 

• Construct Gn from G„+i by contracting the edge X^ and update the 
effective resistances {i?(e)}. 

• Decrease n by one and iterate. 

For sampling uniform spanning trees, more efficient algorithms are known, 
but for the general case the above procedure is the most efficient we are aware 
of. 

We shall need the following lemmas. 

Lemma 21. Suppose {'pk}k=i *5 orthonormal set in L^(A). Then there exists 
a determinantal process with kernel K{x,y) = 4>kix)4>kiy)- 

Proof. For any xi,...,Xn we have {K{xi,Xj))^^^ .^^ = A A*, where Ai^k = 
(j)k{xi). Therefore, det {K{xi,Xj)) is non-negative. Moreover, 


/ det{K{xi,Xj))..dxi...dxn 

J 

= / det {cj)j{xi)). . det {(j).,{xj)).^^dxi.. .dx„ 

J 


sgn(7r)sgn(r) (j)^i^u){xk)(l)r(k){xk) dxi . ..dxn- 


TT^tGSti 


k^l 


In the sum, if7r(fc) ^ T(fc), then J^(j)j^(^i.'j{xk)4>T(k)ixk)dxk = 0, and when 7r(fc) = 
r(fc), this integral is 1. Thus, only the terms with tt = r contribute. We get 



det 


dxi... dxn = ril 
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which along with the non-negativity of det shows that 

^ det is a probability density on A”. If we look at the re¬ 

sulting random variable as a set of unlabeled points in A, we get the desired 
n-point joint intensity pn- 

Lower joint intensities are got by integrating over some of the a;iS: 

Pk{xi,...,Xk) = -^-TTT / Pn{xi,...,Xn)\\.dp{Xj). (13) 

(n — k)\ J ^n,-k ^^j>k 

We caution that m is valid only for a point process that has n points almost 
surely. In general, there is no way to get lower joint intensities from higher ones. 

We now show how to get Pn-i- The others can be found exactly the same 
way, or inductively. Set k = n — 1 in m and expand Pn{xi, ..., Xn) = 
det {K(xi, Xj))^^- as we did before to get 

Pn—l (^ 1 ; ■ ■ ■ 5 Xn—l') 

= '^Sgn{TT)sgn{T)Y[^^^(j}^(^k){xk)^r{k){Xk) / (j)^{n){Xn)^r{n){Xn)dp.{Xn). 

TI-.T d A 

If 7r(n) ^ T(n), the integral vanishes. And if 7r(n) = T(n) = j, tt and r map 
{!,...,n — 1} to {l,2,...,n} — {j} (and the product of the signs of these 
“permutations” is the same as sgn(7r)sgn(r), because 7r(n) = r(n)). This gives 
us 


Pn— 1 (^15 ■ ■ ■j Xn—l ) — det ((/>k(xi))i<i<n-i,k^j det (<(>fe(a:i))fe^j-■ 

We must show that this quantity is equal to det {K{xi, Xj))- j<n-i - ^^is note 
that 

and apply the Cauchy-Binet formula. Recall that for matrices A, B of orders 
m X n and n x m respectively, the Cauchy-Binet formula says 

det(AR) = ^ det {A[ii,..., ini]) det {B[ii,...,im]), (14) 

l<2l ,...,im^n 

where we abuse notation and let A[ii,... ,im] stand for the matrix formed by 
taking the columns numbered ... ,im and B[ii,... im] for the matrix formed 
by the corresponding rows of B. The identity 03 can be proved by observing 
that both sides are multi-linear in the rows of A and in the columns of R. □ 

We now prove Theorem[21 Before that we remark that in many examples the 
kernel K defines a projection operator, i.e, Afc = 1 for all k. Then /fe = 1 for all 
fc, almost surely, and the theorem is trivial. The theorem is applicable to the 
restriction of the process X to D for any Borel set D C A, as the restricted 
process is determinantal with kernel the restriction of AT to D x D. 
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Proof of Theorem^ First assume that /C is a finite-dimensional operator-i.e., 
K{x, y) = ^k4‘k{x)(j)^{y) for some finite n. We show that the processes on 

the left and right side of 0 have the same joint intensities. By m, this implies 
that these processes have the same distribution. 

Note that the process Xj exists by Lemma EH For m > n, the m-point 
joint intensities of both X and Xg are clearly zero. Now consider m < n and 
Xi,..., Xm € A. We claim that: 


E 




det 


To prove itTBl) . note that 


( 15 ) 


{Ki{xi,xj))^^..^^ =AB, ( 16 ) 

where A is the m x n matrix with Ai^p = Ik4’k{xi) and B is the n x m matrix 
with Bk,j = 

Apply Cauchy-Binet formula (d to A, B dehned above and take expecta¬ 
tions. Observe that B[ii,... im] is nonrandom and 

E [det {A[ii ,... ,Zm])] = det {C[ii,.. .,im]) 

where C is the m x n matrix Ci^p = ^k4>k(xi). Applying the Cauchy-Binet 
formula in the reverse direction to C and B, we obtain d and hence also O- 
By Lemma nn given {Ik}k>i, Xg has '^^Ik points, almost surely. Therefore, 

X{A)^Y.^k. 

k 


So far we assumed that the operator K, determined by the kernel K is finite 
dimensional. Now suppose /C is a general trace class operator. Then J2k Afe < oo 
and hence, almost surely, J2k ^k < oo. Therefore, the process Xg is well defined 
and d is valid by the same reasoning. Taking expectations and observing that 
the summands in the Cauchy-Binet formula are non-negative, we obtain 


E 


det 


^ ^ det ... ,2m]) det (.S,..., im ]), 

l<2l 


where C is the same as before. To conclude that the right hand side is just 
det {K{xi, j<m^ Eipply the Cauchy-Binet formula to the finite ap¬ 

proximation {KN{xi,Xj))^^^ j^^, where KN{x,y) = Ylk=i ^k4>k{x)(l)k{y)- Let¬ 
ting N ^ oo, we see that 


E 


det {Kg{x^,xj))^^^^^^^ 


det (iL(a:,,a:j))^<,_^.<^, 


as was required to show. (In short, the proof for the the infinite case is exactly 
the same as before, only we cautiously avoided applying Cauchy-Binet formula 
to the product of two infinite rectangular matrices). 


□ 
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Now we give a probabilistic proof of the following criterion for a Hermitian 
integral kernel to define a determinantal process. 

Theorem 22 (Macchi m , Soshnikov mi). Let K determine a self-adjoint 
integral operator JC on Lp'{K) that is locally trace class. Then K defines a deter¬ 
minantal process on A if and only if all the eigenvalues of IC are in [0,1]. 

Proof. We can assume that K. is trace class, since it suffices to construct a 
determinantal process on compact subsets of A with kernel the restriction of K. 

Sufficiency: If /C is a projection operator, this is precisely LemmaIf the 
eigenvalues are {Afe}, then as in the proof of Theorerndwe construct the process 
Xj. The proof there shows that Xp is determinantal with kernel K. 

Necessity: Suppose that X is determinantal with kernel K. Since the joint 
intensities of X are non-negative, K must be non-negative definite. Now suppose 
that the largest eigenvalue of /C is A > 1. Let Xi be the process obtained by first 
sampling X and then independently deleting each point of X with probability 
1— j. Computing the joint intensities shows that Xi is determinantal with kernel 

Now X has finitely many points (we assumed that K, is trace class) and 
A > 1. Hence, P [A’i(A) = 0] > 0. However, jK has all eigenvalues in [0,1], with 
at least one eigenvalue equal to 1, whence by Theorem 13 ^"[A’i(A) > 1 ] = 1 , a 
contradiction. □ 


Example 23 (Non-measurability of the Bernoullis). A natural question 
that arises from Theorem0is whether, given a realization of the determinantal 
process X, we can determine the values of the Ik’s. This is not always possible, 
i.e., the Ik’s are not measurable w.r.t. the process X in general. 

Consider the graph G with vertices {a,b,c,d} and edges ei = (a, 6 ), 62 = 
{b,c), 63 = {c,d), 64 = {d,a),e 5 = (a, c). By the Burton-Pemantle Theorem @], 
the edge-set of a uniformly chosen spanning tree of G is a determinantal process. 
In this case, the kernel restricted to the set D = {ei, 62 , 63 } is easily computed 
to be 




5 -3 -1 \ 

-3 5 -1 

-1 -1 -ij 


This matrix has eigenvalues 1, ■ Since all measurable events have 

probabilities that are multiples of it follows that the Bernoullis cannot be 
measurable. 


Theorem [7| gives us the distribution of the number of points X (D) in any 
subset of A. Given several regions Dp ,..., Dr, can we find the joint distribution 
of X{Di),..., X{Dr)l It seems that a simple probabilistic description of the 
joint distribution exists only when the Dfs are simultaneously observable, as in 
Theorem El 

Proof of Theorem [3 At first we make the following assumptions: 


• UiDi = A. 
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• K defines a finite dimensional projection operator-i.e., K(x^y) = 
4‘k(x)4>j,(ij) for x,y & K, and {4>k\ is an orthonormal set in L‘^{K). 


Note that by our assumption, (pk are also orthogonal on Di for every 1 < z < r. 
Moreover, it is clear that \(l>k\^■ 

We write 


K{xi,xi) ... K{xi,x„) 

^l) ■ ■ ■ ^n) 

plixi) ... pniXi) \ ( 4>i{xi) ... (piiXn) 
<Pl{Xn) ■■■ 4'n{Xn) j \ (j)„{xi) ... <p.y^{Xn) 


In particular. 


• (17) 


det(is:(a;„a;j))^<._^.<„ = ^ sgn(cr) ^ ^ sgn(r) ^ (a:i) 


^O-GSn 


i=l 


\t^S„ 


i=l 


Now if ki are non-negative integers with ki = n, note that 


(18) 


{X{Di) > ki for all 1 < z < r} = {X{Di) = ki for all 1 < z < r}, 

since by Lemma PTI a determinantal process whose kernel defines a rank-zz 
projection operator has exactly n points, almost surely. Thus, we have 


P [X(Di) = ki for all 1 < z < r] = E 


n 




\ ki 
. 1—1 ^ 


~r~Tl Lr , det{K{xk,xe))T^<f,j,<„dxi...dxr, 

1- • • Kr' D ' ~ ~ 


1 


fci! ■■■krl 


'n: 


D 


n 

^ sgn(cT)sgn(r) ^ (j)^.^{x^)P,Jx^) dxi...dXn. 

(T,r m—1 


Any term with a ^ t vanishes upon integrating. Indeed, if a{rn) ^ T(m) for 
some TO, then 

4>cr^ {Xm)4>T^ (Xm) dXm = 0 

3(m) 

where j (to) is the index for which 



Therefore, 


ki + . . . + < TO < fci -|- . . . -k kjijy^y 




\(l)^^{x)\'^dx 


1 


ki\---kr\ 


n 

(7 m—1 
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Now consider 0 and set M, = for 1 < i < r. Our goal is to compute 

P [Ml = ki,, Mr = kr]. This problem is the same as putting n ball into r 
cells, where the probability for the ball to fall in cell i is Xj^i. To have ki 
balls in cell i for each i, we first take a permutation ct of {1, 2 ,..., n} and then 
put the ball into cell j{m) if fci + ... + < m < ki + ... + kj^jn)- 

However, this counts each assignment of balls 01=1 times. This implies that 


P [Ml = fci,..., Mr = kr\ 



n 

tj m—1 


Thus, 


{XiDi),...,X{Dr)) = iMi,...,Mr), 


(19) 


which is precisely what we wanted to show. 

Now we deal with the two assumptions that we made at the beginning. If 
UiDi 7 ^ A, we could restrict the point process to UiDi. We still have a determi- 
nantal process. Then, if the kernel does not define a projection, apply Theorem[3 
to write T as a mixture of determinantal projection processes. Applying 03 
to each component in the mixture we obtain the theorem. The finite rank as¬ 
sumption can be relaxed in the same way as in Theorem [3 □ 


Applications 

As an application of Theorem[3we can derive the following central limit theorem 
for determinantal processes due to Costin and Lebowitz |3 in case of the sine 
kernel, and due to Soshnikov m for general determinantal processes. 

Theorem 24. Let X„ be a sequence of determinantal processes on A with 
kernels Kn. Let Dn be a sequence of measurable subsets of A such that 
Var {Xn{Dn)) oo as n ^ oo. Then 

X^{ D ) - E [T„ _(j^ ^ ^ ^ 

^Var{Xr,{Dr,)) 

Proof. By Theorem [ 7 | X{Dn) has the same distribution as a sum of indepen¬ 
dent Bernoullis with parameters being the eigenvalues of the integral operators 
associated with Kn restricted to D. A straightforward application of Lindeberg- 
Feller CLT for triangular arrays gives the result. 

□ 

Remark 25. Existing proofs of results of the kind of Theorem |31 (13, |2X|1 
use the moment generating function for particle counts. Indeed, one standard 
way to prove central limit theorems (including the Lindeberg-Feller theorem) 
uses generating functions. The advantage of our proof is that the reason for 
the validity of the CLT is more transparent and a repetition of well known 
computations are avoided. Moreover, by applying the classical theory of sums 
of independent variables, local limit theorems, large deviation principles and 
extreme value asymptotics follow without any extra effort. 
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Radially symmetric processes on the complex plane 

Theorem El implies that when a determinantal process with kernel K has the 
form K{z,w) = 'Yhn^nizw)'^, with respect to a radially symmetric measure /i, 
then the absolute values of the points are independent. More precisely, we have 

Theorem 26. Let X be a determinantal process with kernel K with respect to 
a radially symmetric measure p on C. Write K{z,w) = , where 

QkZ^, 0 < k < n — 1 are the normalized eigenfunctions for K. The following 
construction describes the distribution o/{|zp : z G X}. 

• Let Z be picked from the probability distribution pL/ pl{A), and let Qo = \Z\‘^. 

• For 1 < k < n — 1 let Qk be an independent size-biased version of Qk-i 

(i.e., Qk has density fkid) = with respect to the law of Qk-i)- 

“fc-i 

• Consider the point process in which each point Qk is included with proba¬ 
bility Xk independently of everything else. 

When pL has density </>(|2|), then Qk has density 

TTalq^<p{y/q). (21) 

Theorem 1261 land its higher dimensional analogues) is the only kind of exam¬ 
ple that we know for interesting simultaneously observable counts. 

Proof. Let v be the measure of the squared modulus of a point picked from pL. 
In particular, if y has density (('(l^l), then we have dv{q) = Trifi-y/q) dq. 

For 1 < i < r, let Di be mutually disjoint open annuli centered at 0 with inner 
and outer radii and Ri respectively. Since the functions z^ are orthogonal 
on any annulus centered at zero, it follows that the Dfs are simultaneously 
observable. To compute the eigenvalues, we integrate these functions against the 
restricted kernel; clearly, all terms but one cancel, and we get that for z G Di 

z^Xk,i = / \ka1,{zw)^w^dpL{w), and so 

J Di 

Xk,i = Xkal / \w\‘^'‘dy{w) 

J Di 

fHi 

= ^kal / q''dv{q). 

Jr1 

As ri,Ri change, the last expression remains proportional to the probability 
that the k times size-biased random variable Qk falls in [rl^Lif). When we set 
{ri, Ri) = (0, oo), the result is Xk because OkW^ has norm 1. Thus the constant 
of proportionality equals Xk. The theorem now follows from Proposition El d 

Example 27 (Ginibre Ensemble revisited). Recall that the n*** Ginibre 
ensemble described in Example [El is the determinantal process G„ on C with 
kernel Kn{z, w) = X]fe=o ^kof.{zw)’^ with respect to the complex Gaussian mea¬ 
sure dpi = where a\ = 1/A:!, and Afe = 1. The modulus-squared of a 
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complex Gaussian is a gamma(l, 1) random variable, and its /c-times size-biased 
version has gamma(fc + 1,1) distribution (see (I21ll i. Theorem I2til immediately 
yields the following. 

Theorem 28 (Kostlan ED- The set of absolute values of the points of G„ 
has the same distribution as where Y, are independent and ~ 

gamma(i, 1). 

All of the above holds for n = oo also, in which case we have a determinantal 
process with kernel with respect to dfi = dz. This case is also of 

interest as Goo is a translation invariant process in the plane. 

Example 29 (Zero set of a Gaussian analytic function). Recall the zero 
set of fi( 2 :) := ^ determinantal process in the disk with the 

Bergman kernel 


K{z, w) 


1 

7r(l — zw)^ 


-| 

-^(fc + l)(^^c)^ 

TT 

k^O 


with respect to Lebesgue measure in the unit disk. Theorem applies with 
al = {k + 1) and Afc = 1 (to make K trace class, we first have to restrict it to 
the disk of radius r < 1 and let r ^ 1). From 1121 II we immediately see that Qk 
has beta(A: + 1,1) distribution. Equivalently, we get the following. 

Theorem 30 (Peres and Virag |23p. The set of absolute values of the points 
in the zero set of ti has the same distribution as {t/^ ■ • ■} where 

Ui are i.i.d. uniform[0,l] random variables. 

We can of course consider the determinantal process with the truncated 
Bergman kernel Kn{z,w) = ^J2k=o(^ + l)( 2 :w)^. The set of absolute values 
of this process has the same distribution as ..., 


3. High powers of complex polynomial processes 

Rains m showed that sufficiently high powers of eigenvalues of a random uni¬ 
tary matrix are independent. 

Theorem 31 (Rains HH). Let {zi, ..., Zn} be the set of eigenvalues of a 
random unitary matrix chosen according to Haar measure on U (n). Then for 
every k> n, {zf,... ,z^} has the same distribution as a set of n points chosen 
independently according to uniform measure on the unit circle in the complex 
plane. 

We point out that this theorem holds whenever the angular distribution of the 
points is a trigonometric polynomial. 

Proposition 32. Let {zi,...,Zn} be distributed on (S'^)®” with density 
P(zi,..., z„, zi,... ,z„) w.r.t. uniform measure on (5'^)®", where P is a poly¬ 
nomial of degree d or less in each variable. Then for every k > d the vector 
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has the distribution of n points chosen independently according to 
uniform measure on S^. 

Proof. Fix k > d and consider any joint moment of (zf, • ■ •, z^), 

' n P{zi,...,Zn,Zi,...,Zn)dX, 

(51)»n V / 

where A denotes the uniform measure on (5'^)®". If ^ £i for some i then the 
integral vanishes. To see this, note that the average of a monomial over 
is either 1 or 0 depending on whether the exponent of every Zi matches that 
oi'Zi. Suppose without loss of generality that mi > £i. Then in each term, we 
have an excess of z^ which cannot be matched by an equal power of zi because 
P has degree less than k as a polynomial in zi. 

We conclude that the joint moments are zero unless mi = li for all i. If 
mi = ii for all i, then the expectation equals 1. Thus, the joint moments of 
{zi ,..., z!f) are the same as those of n i.i.d. points chosen uniformly on the unit 
circle. This proves the proposition. □ 

More generally, by conditioning on the absolute values we get the following. 

Corollary 33. Let (^i,... ,fk be complex random variables with distribution 
given by 

P{zi, ...,Zn,Zi,... ,Zn)dlli(zi) ■ ■ ■ 

where P is a polynomial of degree d or less in each variable, and the measures pt 
are radially symmetric. Then for every k > d, the angles Arg(Ci),..., Arg(Cn) 
are independent, have uniform distribution, and are independent of the moduli 
{|Cl|,..-,ICn|}. 

Corollary I2S1 applies to powers of points of determinantal processes with ker¬ 
nels of the form K{z, w) = X]fe=o ^k{zvj)^ w.r.t a radially symmetric measure p 
on the complex plane. Combining this observation with our earlier results on the 
independence of the absolute values of the points, we get the following result. 

Theorem 34. Let X = {zi,..., Zn} be a determinantal process on the complex 
plane with kernel K[z,w) = ^k{zwY (satisfying # {fc > 0 : 0} = nj 

with respect to a radially symmetric measure p. Then for every i > d, the 
points {zf ,..., are distributed as a set of independent random variables with 
radially symmetric distribution. 

4. Permanental processes 

In this section we prove analogous theorems for permanental processes. We begin 
with the following known representation of permanental processes, that can be 
found in Macchi m 

Proposition 35. Let F be a complex Gaussian process on A. Given F, let X 
be a Poisson process in A with intensity jAp. Then X is a permanental process 
with kernel K{x,y) — E [F(a;)F(j/)] . 
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Remark 36. Since any non-negative definite Hermitian kernel is the covariance 
kernel of a complex Gaussian process, it follows that all permanental processes 
are of the above form. 


Proof. Given F, the joint intensities of X are pk{xi,... ,Xk) = 0^=1 
Hence it follows that the unconditional joint intensities of X are pk{xi,... ,Xk) = 

E lull \Fix.r] . Now apply the classical Wick formula (see or m. §2, 
Fact 10) to get the result. □ 


Corollary 37. If K determines a self-adjoint non-negative definite locally trace- 
class integral operator 1C, then there exists a permanental process with kernel K. 

Now we prove Theorem EH using the representation in Proposition 123 We 
need the following simple fact: 


Fact 38. Let 3^ be a Poisson process on A with intensity measure v. Assume 
that v{h.) < oo and ic is absolutely continuous with respect to p,. Let Y be 
the random vector of length 3^(A) obtained from a uniform random ordering 
of the points of y. For A: > 1, the law of Y on the event that 3^(A) = fc is a 
subprobability measure on A^ with density 


gk{zi,...,Zk) 


1 


=-^(A) 


dv 




( 22 ) 


with respect to /r^. We have f gk dp^ = P [3^(A) = k]. 

Proof of Theorem. M (k . We use the construction in Proposition |221 with F{z) = 
where ap are independent standard complex Gaussian ran¬ 
dom variables. Let X be the random vector obtained from a uniform random 
ordering of the points of X. If we first condition on F, then by Fact |2H|the joint 
density of the random vector X on the event {A’(A) = fc} is given by 


jF,k{Zi,...,Zk) 


1 


„-J ,^\Pi^)\^dtJ,(x) 




which is a subprobability measure with total weight P \X (A) = fc | F] . Integrat¬ 
ing over the distribution of F we get that on the event {X{K) = fc} the density 
of X is 


jk{zi,...,Zk) = -^E 


k 




n f ^ ^ \/ m4^m 

i—l \^—^m 


(23) 


which is also a subprobability measure with total weight P [X (A) = fc]. We now 
expand the product inside the expectation 1221) as a sum indexed by ordered 
set partitions (^i,..., Sn) and (Ti,..., r„) of {1,2,..., fc}. The set partitions 
corresponding to a summand q are constructed by letting Si be the set of indices 
i for which q contains the term ^/Xiaicfilzi) and Ti be the set of indices i for 













Hough, Krishnapur, Peres and Virdg/Determinantal Processes and Independence 225 


which q contains the term '/X£a£4>i{zi). The summand corresponding to the 
partitions (Se), (T^) is thus: 


E 


X (n,n 


(n,n 




af ai 


ieSi 




which clearly vanishes unless 15^1 = |T^| for every i. Also note that for a standard 
complex normal random variable a, 


E 


e-^l“l'|o|2™ 


ml 

(1 + A)'"+i ■ 


Therefore by fixing an ordered partition of the integer k with n parts (some 
of the parts may be zero) and then summing over all ordered set partitions 
(Se), (Tf) with those sizes, we find that 


jk{zi,...,Zk) = 77 


s nl, 

rrii—k 

E n,>.n 


A" 


-mi'. 


=i(l + A,)™^+i 


(S.):|S.|=m. 




Now it is easy to see that 





{(t>n{zi) 


(f)l{zk))mi 

4'n{zk)}mn 


(24) 


Recall that the geometric random variables Xi in th® statement of the theorem 
have the distributions P [xi = m]= ■ Therefore we obtain 


jk{zi,...,Zk) 



{^i(2i) 


n " P [Xi = m^] 

(j)i{zk)}mi 


X 


per 


2 


(25) 


\_ {4>n{zi) ••• 4’n{Zk)}mn \ 

Now we integrate (1^ over all the variables Zi. Write the absolute square of the 
permanent on the right as 



{<(> 1 ( 2 : 1 ) • 

4>i{zk)}mi 


{<(> 1 ( 2 : 1 ) • 

(j}i{zk)}mi 

per 

. {<(>n(2:i) • 

■ (t>n{zk)}mn 

per 

. Wn{zi) ■ 

■ 4>nizk)}mn _ 
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and expand these two factors over permuatations tt and a of {1,2,..., A:}. Let- 
ting Ij denote the interval of integers {1 + X]r=i^»’, • ■ ■, J2l=i^r} we get a sum 
of terms of the form 



By orthogonality of (pjs, this term vanishes upon integration unless TT~^{Ij) = 
a~^(Ij) for every 1 < j < n. For a given tt, there are choices of a that 

satisfy this. For each such cr, we get 1 upon integration over ZiS. Summing over 
all kl choices for tt, we get 



p[x{A)=k]= Y. 

rrii—k 
n 

.i^l 


which proves the claim about the number of points in A. Thus by l(23) is a 
mixture of the processes with weights given by nr=iP IXi — ™i]> where 

m = (toi, ..., mn) with rrii being non-negative integers. This is what we wanted 
to prove. □ 


Now we prove Theorem El As before, we remark that it is applicable to the 
restriction of X to D for any Borel set D C A. 

Proof of Theorem MiA Suppose Di, , Dr are simultaneously observable as in 
the statement of the theorem. Use Proposition OH to write X as a Poisson 
process with intensity |U(a:)p where U is a Gaussian process with covariance 
kernel K. Explicitly, F{x) = J2k^k'/^(t>k{x) for x G where ak are 

i.i.d. standard complex Gaussians, i.e., the real and imaginary parts of ak are 
i.i.d. iV(0, i). Then given {ofe}, we know that X{Di), 1 < i < r are indepen¬ 
dent Poisson(/^jF(x)pd^(a:)). Writing JjjjF{x)\^dfi{x) = J2k ^k,i\ak\‘^, we 
see that conditionally given {cfc}, the variables X{Di) for 1 < j < r have the 
same distribution as J2ki'dk,ii ■ ■ ■ iVk,r), where {rjk.ifi^^^r chosen by sam¬ 
pling r]k according to Poisson(Afe|afcp) distribution and then assigning rjk to the 
cells Di multinomially with probabilities 

Now when we integrate over the randomness in {ofe}, we see that r]k has a 
Geometric distribution with mean Afc and given rik, the vector {rik,i, ■. ■, rjk.r) is 
still Multinomial ■ ■ ■, This completes the proof. □ 


Generalization: a-determinantal processes 

One way to generalize the concept of determinantal and permanental processes 
is to consider point processes with joint intensities given by 

n 

Pn{xi,... ,Xn) = deta{K{xi,Xj))‘^= ^ ,a;,r(i)), (26) 
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where v{t^) is the number of cycles in the permutation tt. 

Such point processes are called a-determinantal processes and were intro¬ 
duced by Vere-Jones in m . The values a = —1 and a = -1-1 correspond to 
determinantal and permanental processes, respectively. It is easy to check that 
the proof of Theorem dean be modified to get: 

Proposition 39. For a point process with joint intensities given by \2f)\) (when 
it exists), the number of points that fall in any subregion is: 

• a sum of independent Binomial(— — ,—aXk) random variables, if—— is a 
positive integer. 

• a sum of independent Negative Binomial(^, random variables, if 

a > 0. 

In fact, if — — is a positive integer, this process is just a union of — — i.i.d. copies 
of the determinantal process with kernel —aK. Similarly, if ^ is a positive in¬ 
teger, this process is a union of ^ i.i.d. copies of the permanental process with 
kernel aK. More generally, the union of m i.i.d. copies of the process corre¬ 
sponding to a and kernel K gives a process distributed according to ^ and 
kernel mK. If K is real, then ^-determinantal processes also exist m For 
a > 0, little is known about the existence of a-determinantal processes beyond 
these examples. Shirai and Takahashi conjecture the following: 

Conjecture 40. If K is a (real) kernel defining a self-adjoint, non-negative 
integral operator on and 0 < a < 2, then the a-determinantal process 

with kernel K exists. However, if a > 2, then there exists such a kernel K for 
which there is no corresponding a-determinantal process. 

We verify this conjecture for a > 4. Let A be a discrete space consisting of 
three points, and consider the 3x3 matrix 

2 -I -1 \ 

-1 2 -1 . 

-1 -I 2 / 

It is easy to check that the eigenvalues of K are 3,3, 0 and 

deta{K{i,j))i<ij <3 = 2(4 - a)(a -b I), 

which is negative for a > 4. Since point processes must have non-negative joint 
intensities, we conclude that no a-determinantal processes with this kernel can 
exist for a > 4. 

Remark 41. It was pointed out to us by Steve Evans that this question is 
closely related to the question of inhnite divisibility of a multivariate version of 
Gamma distribution that has been studied by Griffiths m and Griffiths and 
Milne m In short, knowing that a certain multivariate Gamma distribution is 
not infinitely divisible would tell us that a-permanental processes for a related 
kernel do not exist for sufficiently large a (but does not give us the smallest a 
for which this happens). 
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